Superconductivity and Abelian Chiral Anomalies 



Y. HatsugaO and S. Ryu 

Department of Applied Physics, University of Tokyo, 
7-3-1, Hongo, Bunkyo-ku, Tokyo 113-8656, JAPAN 

M. Kohmoto 

Institute for Solid State Physics, University of Tokyo, 
5-1-5, Kashiwanoha, Kashiwa, Chiba, 277-8581, JAPAN 
(Dated: Nov. 11, 2003) 

Motivated by the geometric character of spin Hall conductance, the topological invariants of 
generic superconductivity are discussed based on the Bogoliuvov-de Gennes equation on lattices. 
They are given by the Chern numbers of degenerate condensate bands for unitary order, which 
are realizations of Abelian chiral anomalies for non-Abelian connections. The three types of Chern 
numbers for the x, y and z-directions are given by covering degrees of some doubled surfaces around 
the Dirac monopoles. For nonunitary states, several topological invariants are defined by analyzing 
the so-called g-helicity. Topological origins of the nodal structures of superconducting gaps are also 
discussed. 



The importance of quantum-mechanical phases in con- 
densed matter physics has been recognized and empha- 
sized for recent several decades. The fundamental char- 
acter of a vector potential is evident in the Aharonov- 
Bohm effect where the C/(l) gauge structure is essential 
and a magnetic field in itself plays only a secondary role 
0. Topological structures in quantum gauge field theo- 
ries have also been studied and extensive knowledge has 
been accumulated 0- Quantum mechanics itself supplies 
a fundamental gauge structure^]. It is known as geo- 
metrical phases in many different contexts, where gauge 
structures emerge by restricting physical spaces. The 
quantum Hall effect is one of the key phenomena to estab- 
lish the importance of geometrical phases The topo- 
logical character of the Hall conductance was first real- 
ized by the Chern number expression, where the Bloch 
functions define "vector potentials" in the magnetic Bril- 
louin zone accompanied with a novel gauge structure!^]. 
Further the ground state of the fractional quantum Hall 
effect is a complex many-body state where another kind 
of gauge structure emerges (6j. These quantum states 
with nontrivial geometrical phases are characterized by 
topological orders which extend an idea of order parame- 
ters in statistical mechanics to the quantum states with- 
out spontaneous symmetry breaking 0. Also we point 
out an importance of boundary effects for topologically 
nontrivial systems. Bulk properties are closely related 
to edge states and localized states near impurities and 
vorticesHliElEilEl. 



Topologically nontrivial structures in superconductors 
also have a long history. Recently, following a prediction 
of flux phases for correlated electron systems ; 13]. spin 
Hall conductance is defined for superconductors based 
on the Bogoliuvov-de Gennes (B-dG) equation 0,0,0] 
[sij l . As for singlet states and triplet states besides equal- 
spin-pairing states, a map to a parameter space which 
represents the BdG hamiltonian is considered. 0| In the 



parameter space, an analogue of the Dirac monopole [l^ 
exists and the Chern numbers are analyzed. 

In this paper, we establish a topological characteri- 
zation of general superconductivity based on the BdG 
equation on lattices. The energy spectrum of the BdG 
hamiltonian are fully used to calculate the Chern num- 
bers of the superconductors. As for the unitary su- 
perconductors, condensed matter realizations of chiral 
anomalies for non-Abelian connections are given explic- 
itly. Topological consideration is useful to distinguish su- 
perconductivities with the same pairing symmetry. The 
present analysis also clarifies nodal structures of super- 
conducting gaps with various anisotropic order param- 
eters, which is closely related to the quantum Hall ef- 
fect in three-dimensions 0, H^]. Various types of the 
nodal structures are not accidental but have fundamental 
topological origins. A possible time-reversal symmetry- 
breaking and an unconventional gap structure are pro- 
posed based on the experiments |33| 

Bogoliuvov-de Gennes hamiltonian. Let us start from 
the following hamiltonian on lattices with spin-rotation 
symmetry: 



where Cia is the electron annihilation operator with 
spin a at site i, = t*,, tA-i-^--3-4 ^ (^^j4-3;-2^i 
ya^a^.a^a^ ^ y^<T.<ri ;<t4^3 ^^^^ ^ -^^ ^ chcmical potential. 
Summations over repeated spin indices a are implied 
hereafter. The mean field BCS approximation leads to. 



^ =^^^j4aCj^ + ^{^1f''*Cja^C^„^ + h.C.) - ^ 



where the order parameters {'^if — —^ji are given 
by A-^p = l^J^"'="'"'(cj.3C,„J. The usual mean field 
theory leads to the gap equation of which a solution 
gives an order parameter. Here we do not follow this 
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procedure but a priori assume order parameters which 
may be reaUzed for some interactions V^'^^'^^'^^'^''. Let 
us consider the two cases separately ^21, 22|: (i) singlet 



states A^, 
states A, 



ipijidy, {tpij = ipyj) and (ii) triplet 



[dij ■ (T)iay, {di. 



-dji), where 



{AijY'^ — A°j^ is a 2 X 2 matrix in the spin space and 
denotes matrix transpose. Now assume the translational 



symmetry, namely, ti 



t{i - j), = - j) and 



also the absence of a magnetic field, that is, t{i — j) to 
be real. Then, except a constant, the BdG hamiltonian 
is given by a 4 x 4 matrix hk as 



H=^clhkCk, hk = 



efeCTo Afc 



where c| = (c|(A;), c|(fc), C|(— fc), Cj(— fc)) with Ca{k) = 



ik-r. 



Afe = 
and V 



E,e-^'^-^^A(£), A_, = -A,, = (^J J 

is a volume of the sysmtem. The order parameter is 
given by A^ = ipkicTy, A^ = — A^ for singlet states and 
Afe = {dk ■(T)i(jy, Afc — Afc for triplet states ( tpk is even 
and dk is odd in k). |^ 

The BdG hamiltonian has a particle-hole symmetry. 

"^■"l = Ekf'"'), then af"^- 

eigenstate with energy —Ek where C — pxK for singlet 
states and C — —ipyK for triplet states ( Uk and Vk are 
the two-component vectors and K is, a, complex conjugate 
operator and the Pauli matrices p operate on the two 
component blocks.) 13 Then it is useful to consider 



If hk 



is also an 



hi 



4po 



AkAl 







AtA, 



For singlet states, we 



have AkAj. = Aj,Ak — IV'fcPo'o and for triplet states. 



^kAl 



k - idkxd*k, 



<ife I o'o + Q/c ■ with a real vector q 
which we call g-helicity ( f represents hermite conjugate 
and * complex conjugate). 

Chern numbers for unitary states. Singlet order and 
triplet order with vanishing g-helicity are called uni- 
tary since A^A^ — Aj^A^ oc ctq. Nonunitary triplet 
states (Qk ^ 0) will be discussed later. For unitary 
states, we define a unitary matrix A° by A^ — |Afc|A°, 
(|Afe| = \ipk\ for singlet states and |Afc| = \dk\ for triplet 
states, respectively). Since the spectra are doubly degen- 
erate as will be shown later, fixing phases of the states 
is not enough to determine Chern numbers by the stan- 
dard procedure 0, 0| . Instead, one can define non- 
Abelian vector potentials and fluxes following definitions 
of generalized non-Abelian connections (2^ . 

Let us assume that the states are M-fold degener- 
ate { M = 2 in the present unitary case) as |a), 
a = 1, ■ • • , M. Then a non-Abelian connection is de- 
fined by Af = (a\d^,\P), ^"'^ = Afdkf, and 9^ = 
(9fe^,/i = a;,?/,z.|33 A unit 



"a gauge transformation" A. = u!^ Au> + uj^ duj . Then the 
field strength — dA. -|- A A ^ is gauge covariant since 
= iJ^Tuj. One may write T — ^Ff^^dkf^Adk,^, F — 
d,,Au - d^A^ + [A^,A,], {{A^f^ = Af ). Then "a 
magnetic field" in the /i-direction is — ^e^jyATr 
Since Tr^ is unitary invariant, so is B^. The total fiux 
passing through the i/A-plane is given by an integral of 
the magnetic field B^ over the two-dimensional Brillouin 
zone, T"^^, ( is fixed). The first Chern number is 

X./^-^=4i l^Jk.^dk.B, 

This is the Abelian chiral anomaly discussed in the non- 
Abelian gauge theories 0, Here we have considered 
the cubic lattice. Extensions to other lattice structures 
is straightforward. 

Dirac monopole in the parameter space. 
The BdG equation for the unitary states 
efefJo |A.|A2 \ ( Uk\ f Uk 



\Ak\iK) 

duces to 



0\-l 



|A.| 



lAfel 



Vk 



Uk 



Ek 



Vk 



=Ek 



Uk 

Kvu 



Thus energies are given by Ek = ±R {R = \/ e^ + |Afc|^) 
and the states are doubly degenerate. The band with 
energy —R is the superconducting condensates of pairs. 
On the other hand, the band with energy +R repre- 
sents quasiparticle excitations. By a parametrization: 
ek — Rcos9 and \Ak\ — Rs'mO, eigenvectors of conden- 

Uk 

I — I ^f r/ )^ 



sate (Ek = —R) are i . n 



a) , where 



\R,o)p = 



cos ■ 



and \a)a-, {a ~ 1,2) are arbitrary 



orthonormalized states in the cr-space. Now let us take 
eigenstates of the A^ with eigenvalues e~^'^" {a = 1, 2) 
to calculate the Chern numbers. 37] The degenerate or- 
thonormal eigenvetor for the condensate band Ek = —R 

The connection is given by 



is given by \^a) = 



e^'t'" cos f I a) 



A'f ^{iPo.\dMp) = {R^\d^,\Rp)MP)- + {Rc^Rp)M^^. 
=Afip)6^0 + {R^\Rp)pAf{a), 

where Af{p) = with \Rc,)p = 

\R,0,^o.)p = ( Z'Zh)- Also Af {a) = (a|a^|/3). = 



generate state \a) — > \a) 



{Wdf,U)ap with U = Then the total mag- 

netic field in the parameter space is B^ = B^{p) + 
B^{a), where B^(p) = e^^ATr i9^Aa(p) and B^(cr) = 
^tit^x^T: duA\{a). Since (a) vanishes by the "sum rule" 
among the filled bands, we have B^ ~ -S^(p). It im- 
[ijj: unitary) causes plies that the Chern numbers, C^^^, of the condensed 



unitary transformation of a de- 
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band in the /Lt-direction, is given by the sum of the Chern 
numbers of the two vectors \Ra)p, {a = 1,2) which are 
the eigenstates of the 2x2 hamiltonians 

where Ra.x = R 

sill ^ COS Ra.Y — R sin SITL (I)q, and 

Ra,z = Rcos9. Namely they are — X^q^^.m- ^ow 
we have reduced the problem to calculate the Chern num- 
bers of the eigenstates of the 2x2 matrices, h'^. By 
mapping from the two-dimensional Brillouin zone to the 
three-dimensional space, T^-^ 9 {k^,k\) — > Ra, we ob- 
tain a closed oriented surface Ra{T^)^). The wrapping 
degree of the map around the origin gives a charge of the 
Dirac monopole sitting there. This is the Chern number 

For the present degenerate case, the map from a two- 
dimensional point to two three-dimensional points T^^ 3 
{kv, k\) {i?Q=i, Ra=2\ defines ( fixing k^) the surfaces 
{Ra=i{T'^x),Ra=2{T'^x)} which determine the two cov- 
ering degrees of the maps around the origins, Na-^^\lk^), 
[a — 1,2). They give the Chern numbers Cq,,^ respec- 
tively. Since only the condensed states are filled for the 
superconducting ground state, the Chern numbers of the 
unitary states are given by 

The Chern numbers defined here for the unitary super- 
conductors satisfy C^\{ki_i) = 4 x integer for the singlet 
order and Cy\{kfj) = 2 x integer for the triplet order. |33 
Nonunitary states. In these triplet states, there are 
no degeneracies in solutions of the BdG equation. There 
are four quasiparticle bands, which are classified by the 
q-helicity as 

(o- • 9fc)M± = ±qkU±, qk = \qk\, 
hl^f^{el + \dkf±qk)^jf, 

"^^^ " ( "± ) ' " ( ) ' ^± " -i<^yU^- 

Then states with helicity +qk and energy ±E+q — 
±\J(\ -t- Idfep + gfe, are \ ±E+q) = U +7]%^, where U+ = 

I I \ orthonormal vectors ri^„ are 

V2 \V+ -V+ J '+« 

determined as the eigenvectors of the reduced 2x2 hamil- 
tonian h+g — U\ghU+q with energies ±-B+g- The 
hamiltonian h+q is traceless and it can be expressed 
by the Pauli matrices, h^q = a ■ R+q where R+q = 

(e,Im(i_| ,Re(i^ ) is a real vector and = M^(<ife • 

a)u-. As for the helicity —qk state, one can follow al- 
most the same procedure and obtain the reduced BdG 
hamiltonian similarly. Here we can define several topo- 
logical invariants. As discussed above, the states with 



zbg^-helicities and energies ±i?±g are non-degenerate, 

where E^^q{k) = i/lefep + |dfc|2 + Sq, {Sq = ±qk and 
Ee = i)- Then the standard Chern number C'^{kp) in 
the /Lt-direction for a fixed fc^ is obtained by the standard 
wayflj. 

Also we have topological invariants in the /Lt-direction 
N'i(k^j). They are wrapping degrees around the ori- 
gin of the map {k^,kx) Q{k) — q^, which define 
closed surfaces q{T^x) three-dimensions [l7| . Further 
we have other topological invariants iV^(/e^), which are 
also wrapping degrees of the map around the origin, 
{ki,,kx) — > R±q{k). (The reduced hamiltonians are 
h±q — R±q ■ cr for the g-helicity ig^.) 

Up to this point, topological arguments have been ap- 
plied to gapful cases. However, the nodal structure of 
the gap function is, in fact, characterized by the topo- 
logical description. Formally we have treated a three- 
dimensional superconductivity as a collection of two- 
dimensional systems parametrized by, say, kz- In R- 
space, closed surfaces parametrized by (kx,ky) are gener- 
ically away from the monopole at the origin. As kz is var- 
ied, the surfaces move around and they can pass through 
the monopole. Since the distance between a point on 
a surface and the monopole gives a half of the energy 
gap Eg{kx, ky]kz), the gap closes at the value of kz when 
the monopole is on the surface. Thus the nodal struc- 
ture of the superconductivity is point-like generically if 
the time-reversal symmetry is broken. When the Chern 
number jumps as kz varies, the superconducting gap has 
to be closed due to topological stability of Chern num- 
bers. Also a non-zero Chern number implies that the 
corresponding two-dimensional system has a non-trivial 
topological order. The superconducting node is consid- 
ered as the critical point of the quantum phase transition 
between two states with different toopological orders. 

To make the discussion clear, let us take an exam- 
ple Ak ~ ci°(sinfca; -I- isinfc,,)cr;r, = — 2i(cosfcz -I- 
cosfcy -I- COS kz) — fi, {t > 0) jrgl- This is an analogue 
of the Anderson-Brinkman-Morel(ABM) state in ^He su- 
perfluid. For a fixed value of fcz, the surface is reduced 
to that of the chiral p-wave order parameter with a mod- 
ified chemical potential, /j, — 2t cos kz . ( we can recover 
the ABM state, dk {0,0, d^ikx + iky)) in the limit 
of /i —6t + 0,k 0.) There are two quantum phase 
transitions which are accompanied by jumps of the Chern 
number between —2 and for —6t < fi < ~2t. These cor- 
respond to the gap nodes at the north and south poles 
on the Fermi surface. In Fig. 1, surfaces Ri{T^y) are 
shown with the monopole at the origin. The Chern num- 
ber jumps between —2 and +2 for — 2t < /.i < 2t and 
and +2 for 2t < /i < 6t. 
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(a) (b) 




FIG. 1: Examples of closed surfaces Ri {T^y) which are cut by 
the _Ry — i?z-plane. The monopole is at the origin, t = = 1, 
11 = -5: (a)fc3 = 0, (b)fc^ = -2tv/5. 

For line nodes, we need additional constraints to keep 
the monopole on the closed surfaces when is varied. 
To make a discussion simple, we take singlet order or 
triplet order with dx — dy = and dz ^ 0. Further let 
us require that the order parameters are real, namely we 
have a time-reversal symmetry. Then the closed surfaces 
in the i?-space collapse into a board like region on the 
Rx-Rz plane and one can expect a situation where the 
monopole moves along the surface when is varied. [25j 
Thus a line node appears in the superconducting gap. As 
shown in this example, the nodal structure of the super- 
conductivity has a fundamental relation to topological 
order. The detailed discussion on this point will be given 
elsewhere. |2 6] 
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